Analytic solution for relativistic transverse 
flow at the softest point 



Tamas S. BiroQ, 

", MTA RMKI, H-1525 Budapest, P.O.Box 49 

PACS: 25.75 



O 

^ ■ Keywords: transverse flow, analytic solution 

^ : 

' Abstract: We obtain an extension of Bjorken's 1+1 dimensional scaling relativistic flow solution 

^ ■ to relativistic transverse velocities with cylindrical symmetry in 1+3 dimension at constant, 

Tlj" i homogeneous pressure (vanishing sound velocity). This can be the situation during a first order 

I phase transition converting quark matter into hadron matter in relativistic heavy ion collisions. 

0\ ■ Hydrodynamics often allow for nonrelativistic scaling solutions. Relativistic flow, however, 

^ seems to be an exception: besides Bjorken's 1+1 dimensional ansatz and the spherically 
symmetric relativistic expansion, no analytical solution is known [^. 

O I In this paper we present an extension of Bjorken's ansatz p| for longitudinally and 



X 



transversally relativistic flow patterns with cylindrical symmetry in 1+3 dimensions. This 
is an analytical solution of the flow equations of a perfect fluid for physical situations when 
the sound velocity is zero, cj. = dp/de = 0, with energy density e and pressure p(e). 

In particular this happens during a flrst order phase transition, the pressure is constant 
while the energy density changes (in heavy ion collisons increases and drops again). This 
should, in principle, be signalled by a vanishing sound velocity. A remnant of this effect 
in flnite size, flnite time transitions might be a softest point of the equation of state, 
where cj. is minimal. In fact, this has been suggested as a signal of phase transition by 
Shuryak0, and investigated numerically in several recent works 0, |^. 

In the light of this research, the presentation of an analytical solution including relativistic 
transverse flow is worthwile. For a general equation of state p{e) the analytical solution 
given in this paper cannot be extended, only a perturbative expansion in terms of mild 
(v -C 1) transverse velocities can be established. Such an approximation has been recently 
presented in [^. Nonrelativistic analytic solution has been also given several times, with 
respect to heavy ions see 0, |^ . 



The 1+1 dimensional Bjorken flow four-velocity is a normalized, timelike vector. It is 
natural to choose this as the flrst of our comoving frame basis vectors (vierbein). The 
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three further, spacehke vectors will be constructed orthogonal to this, separating the 
two transverse directions. This basis fits excellently to a cylindrical symmetry and to 
longitudinally extreme relativistic fiow. 

<i={o.oAK}, e» = (o,o,-^,i), (1) 

with t time coordinate and z longitudinal (beam-along) coordinate, x and y transverse, 
cartesian coordinates. The cylindrical radius is given by, r = \^x^ + y^, and r is the 
longitudinal proper time: r = y/t'^ — z^. 

This way the basis can be re-written in terms of a hyperbolic angle (coordinate-rapidity) 
77 and a polar angle 

e° = (cosh T]^ sinh 77, 0, 0) , ~ (sinh 77, cosh 77, 0, 0) , 

= (0, 0, COS0, sin 0, ) , = (0, 0, — sin 0, cos0, ) . (2) 

This basis is orthonormal, • e^'^ — g°^^ and satisfies the differential relations: 

del = eldr], de^ = ejd7;, del = eld(j), del = -eld(j). (3) 

The space-time coordinate differential and the partial derivative are in our basis given by, 

dXf^ — e^^dr + e^^rdr] + e^dr + e^^rdcf), 

d = e'^±-e'---e'--e'-— U) 
^ ''dr ^rdr] "dr ^rd(t>' ^' 



We consider ideal ffuids ("dry water"), where the energy momentum tensor is given by 

T^,u = (e + p)u^,u^ - pg^,u, (5) 

and the equation of state is given in the form of p(e). The ansatz for an almost boost 
invariant ffow with some transverse, cylindrically symmetric component is then given by 

Un^l (ej + vefj , (6) 

using the Lorentz factor 7 = (1 — v^)"^^'^. This four- velocity is normalized to one: 

u^u^" = 7^ - 7^^^ = 1. (7) 

For the four divergence of the ffow we obtain 

5X-^7 + ^(7t')+7(- + -j. (8) 
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The co-moving flow derivative is 

In order to formulate the Euler equation describing a possible acceleretion or decceleration 
of the flow, one also needs a projection orthogonal to u^. We use 

Va = ((7am - uxu^) d", (10) 

which is given by 




The relativistic equation of local energy and momentum conservation, d^T^^ = 0, can 
then be projected to a component parallel to and components orthogonal to that. 
Denoting the co-moving derivative (H) by 7 times an overdot, (i.e. u^d^f = 7/ for any 
/), we arrive at 

7e + wdf,u^ = 0, 



w-fUx-VxP = 0. (12) 



Here we introduced the enthalpy density w = e + p. The first equation of (|T^) is the 
local form of the dE + pdV = adiabatic flow condition; in fact this is equivalent to 
the conservation of the entropy flow in a one-component matter. The second equation 
resembles the familiar form of the nonrelativistic Euler equation. 

Summarizing we have alltogether four independent equations, 

\OT or J 
d d \ ( d d 



\OT or J \0r or J 

1 d 

-^P = 0, 



T drj 

Id 
r 0(p 

Here the four divergence (^ without the Lorentz factor 7 is denoted hj D = ^d^u^, which 
can be expanded to 

2 f 9 d \ 1 V 

^ = ^ hr + ^7r P + - + -- 14 

\or OT J T r 

Let us restrict ourselves now to p = po = constant situations. In this case = and 
w = e+po- From (|l^) we arrive at two independent equations only, the first and the second 
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one. Since derivatives of the pressure vanish, we obtain an equation for the transverse 
flow component alone: 

(^ + >'|)-0. (15) 

This equation is valid even for relativistic transverse flow, it is a quasilinear partial differ- 
ential equation. A general solution can be obtained by the factorizing ansatz v = a{T)-b{r). 
The above eq.(|^) leads to 

a'ir) 



aHr) 



+ b'{r) = 0. (16) 



Here both terms are constant, balancing each other to zero. The solution, which is regular 
at the cylindrical axis r = is given by 



ar 



1 + ar 

The first equation of the system (|T3D is the coohng equation. It reduces to 

f d d\ 



(17) 



with D given by eq.(n). Utilizing the solution (K% for v{T,r) we obtain 



D = i + -^, (19) 

r 1 + ar 

which can also be written as a comoving derivative, 
with 

E(r) = r(l + ar)l (21) 
Using D in this form the cooling equation (|1^) can be rewritten as 

+ + = ^^^^ 

A particular solution of this equation is, when the quantity on which + v-^^ operates 
is constant, 

const. 
r(l + ar)^ 

This solution interpolates between the one- dimensional cooling law e oc 1/r and a three- 
dimensional one for large longitudinal proper time, e oc 1/r^. 

There is, however, a more general solution to the "cooling law". First inserting the 
analytic solution (0) for v in the co-moving derivative operation we get, 

(i- + T^^l"! + log^(^)) = 0- (24) 

\oT 1 + ar or I 



Multiplying this by (1 + ar)/a we obtain 



+ (logw + logE(r)) = 0. (25) 



\d\og{l + ar) d\ogr 
This equation is of type 

whose solution is simply a function of a; — y, -F(a;, ?/)soiution = F{x — y). Henceforth we 
obtain 

logw + logE(r) = F(logr-log(l + ar)). (27) 

In this result one easily recognizes the analytic form of the transverse flow velocity v ( [T7| ) 
after collecting the terms in the argument of the general function F inside one logarithm. 
Finally we arrive at 

\ogw + log E{t) = F{v/ a) = log f{v) (28) 
as an analytic solution of the cooling law with relativistic cylindrical transverse flow: 



e[r, T) = -po + 



r(l + ar 

a 



1 + ar 



r. (29) 



The unknown function of one variable f{v) can be used to match the initial profile of the 
radial distribution of the energy density at r = tq, 

e{r, To) = -Po + — -J (-r^) ■ (30) 



Let us finally discuss some properties of this analytic solution. First one would like to be 
convinced that the Bjorken scaling limit is contained in eq.(p9|). (The 1+1 dimensional 
solution was presented by Gyulassy and Matsui in 1984 [0.) This is indeed the case, for 
f — > one obtains a = and 

e(r,r) = -po + ^, (31) 
r 

leading to 

(e(r)+po) = (eo+Po)-. (32) 

r 

As an example of the p = const, situation let us consider an oversimplified system: mass- 
less quark-gluon plasma described by the bag equation of state, 

eg = aT^ + B, pq = ^aT^ - B, (33) 

keeps a Gibbs equilibrium with a light, relativistic pion gas, described by 

en = hT\ pn = ^hT\ (34) 
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During a first order phase transition, for simplicity assumed to take place in the total 
volume, the pressure is constant and homogeneous (Gibbs criterion). Actually from this 
requirement pq = pn one usually obtains the transition temperature due to 

^-ia-h)T^ = B. (35) 

The energy density is that of a mixture, containing x part quark matter and (1 — x) part 
hadron matter, 

e = xeQ{T,) + {l-x)eH{T,). (36) 

This partition x is what changes according to the cooling law, respectively its solution 
(pQl). Utilizing eqs.( p3| , p4D and (^) we arrive at 

e = 4Bx{t) + eH = -pn + (cq + Pg) (37) 

r 

Assuming at r = tq pure quark matter x{tq) = 1, the time of the total conversion ri when 
x(ri) = is given by 

Ti = ro = -To. (38) 

It is determined by the relative number of degrees of freedom in the two phases in this 
simple scenario. 

The same discussion is somewhat more complex using the analytic solution with a cylin- 
drical flow. At the beginning of the phase transition x{tq) = 1 let be a transverse velocity, 

r 

v{To,r) = vo — . (39) 
/to 

The analytic solution (^) leads to 

v(T,r) = — ^ -r (40) 

Ro + Voir-To) 

featuring a slowing down of the transverse flow. From the initial condition for the trans- 
verse flow we also obtained the parameter a 

Since the four-flow is given by 

«M = 7 (ej + vel) 

a slowing down of the transverse component implies a slowing down of the original Bjorken 
component as well (7 is decreasing to one, when v is decreasing to 0). 

The cooling of the energy density can be obtained now by fitting the profile at r = tq to 
e(r,ro) = e(r), 

e(r,r)+po = ^§(e(f r)+po) (42) 

with 

RiT)=Ro + Vo{T-To). (43) 
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This solution follows the Bjorken coohng law at the beginning r tq, 

e(r, r) + po ~ — (e(r) + po) , (44) 

T 

but for long times r ^ tq turns over to a three-dimensional scaling 

e('-.-)+P.-^f (e(^^)+P.). (45) 

Here for large times, the driving force besides the constant pressure will be e(0), the 
original energy density at the axis. 

Finally wc repeat the calculation in the simplified scenario assuming a first order phase 
transition between ideal quark gluon plasma and ideal pion gas. The general solution for 
the quark matter part, x is given by 

Assuming a linearly decreasing transverse profile of the QGP part initially (at the begin- 
ning of the phase transition) , 

x{To,r)^l-r/R,, (47) 
we obtain the transverse radius of the mixed phase as 

r,^^RU-J^(l^-l]). (48) 
Ro \ a-h\ToRl J J ^ ' 

In the simple case, when Rq = vqTo, the scaling transverse radius R is proportional to 
the longitudinal proper time t: R — vqt. The transverse velocity field is also particularly 
simple: v — r/r. The above result simplifies to 

n^Ri-[ r r- • 49 

To \a — h a — Htq J 

This expression initially grows, achieves a maximum and then decreases towards zero. 
The maximum, 

is achieved at time 

a \i/3 



^^^^^'^AhJ ' ^^^^ 



the a; = for the whole space (equivalent to Vb = 0) at 

Realistic estimates use a = 37 for the quark gluon plasma and h = 3 for the pion gas. 
This leads to r™'^'' ^ 1.14i?i, r™'^'' ^ 1.45ro. and Ti ^ 2.31ro. It is realistic to assume 
To ~ 5 fm/c for a CERN SPS Pb-|-Pb experiment. The radial extension of the quark 
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matter can grow with about 14 per cent and then the conversion into hadrons eats it up, 
reaching the zero radius in about 11.5 fm/c. This is a fast hadronization even in this 
simple scenario. 

For a smaller initial radial flow, realistically Vq = 0.6 at the radial edge of the cylinder 
r = -Ro) one obtains VqTq/Rq = 0.428. Fig.l shows the radial space - time evolution 
profiles of the mixture: the outmost curve corresponds to x = (pure hadron matter) 
and the steps are 0.1. The innermost curve belonging to x = 1.0 shrinks to a point, since 
the phase conversion starts immediately at r = tq. Here we assumed a linear initial profile 
x(r,ro) = 1 - r/Ro. 

With these parameters the mixed phase of Rq = 7 fm longs for about ~ 17.5 fm/c. 

In conclusion we have presented an analytical solution for relativistic transverse flow of a 
perfect fluid at the softest point. This can be reahzed during a first order phase transition, 
when the pressure is homogeneous and constant in a large volume. The phase conversion 
during this stage of expansion is given by this analytical solution which scales with the 
longitudinal proper time initially like the Bjorken flow 1/r, but eventually like a spherical 
flow 1/r^. For simple equations of state for the quark and hadronic matter side the 
estimated longitudinal proper time spent in the mixed phase turns out to be more than 
twice the time when the phase transition began. 
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Figure 1: Contour lines of several QGP - hadronic mixtures (from to 1 by 0.1 steps) in 
the t/tq - r/Ro plane. The initial profile at r = tq was assumed to be linear giving zero 
percent quark matter exactly at Ri = Rq. The initial transverse flow is Vq = 0.6 at the 
edge R = Rq. 
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